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Abstract 


Ground  motion  records  from  seven  high  explosive  cratering  events 
in  northeastern  Montana  were  analyzed  for  peak  velocity,  power  spectral 
density,  and  velocity  spectra.  The  events  included  four  20 -ton  single 
charges  at  depths  of  burst  which  varied  from  42  to  57  ft,  a  140 -ton  row 
charge  consisting  of  three  20 -ton  charges  and  two  40 -ton  charges  at 
optimum  depths  of  burst,  and  a  fully  coupled  charge  of  0.5  tons  and  a  de¬ 
coupled  charge  of  0.5  tons  at  optimum  depths  of  burst.  It  was  found  that 
at  these  depths  and  charge  weights  an  increase  in  depth  of  burst  resulted 
in  an  increase  in  peak  velocities  and  power  spectral  densities  as  measured 
at  distant  points  (>5  km),  while  no  significant  frequency  shifts  were  noted. 
Power  spectral  density  was  found  to  be  approximately  proportional  to  the 
first  power  of  yield.  For  this  region  it  was  determined  that  power  spectral 
densities  varied  inversely  as  radius  to  the  3.55  power,  and  peak  velocities 
varied  inversely  as  radius  to  the  1.6  power.  An  increase  in  both  velocities 
and  power  spectral  densities  for  small  decoupling  factors  was  found  to 
occur  for  a  certain  explosive -cavity  configuration.  Three  analysis  tech¬ 
niques,  peak  velocity,  velocity  spectra,  and  power  spectral  density,  are 
compared  and  it  is  shown  that  power  spectral  density  is  the  most  consistent 
method  when  comparing  records  from  different  measuring  stations. 
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INTERMEDIATE  RANGE  GROUND  MOTIONS 
FOR  PRE-GONDOLA  II  AND 
ASSOCIATED  EVENTS 

Introduction 


DESCRIPTION  OF  PROJECT  PRE- 
GONDOLA  II 

Project  Pre -Gondola  II  was  a  row- 
charge  cratering  experiment  in  weak,  wet 
clay-shale  conducted  by  the  U.  S.  Army 
Engineer  Nuclear  Cratering  Group  (NCG) 
as  a  part  of  the  joint  Atomic  Energy 
Commission -Corps  of  Engineers  nuclear 
excavation  research  program.  The  pri¬ 
mary  purpose  of  this  nominal  140 -ton 
five -charge  experiment  was  to  gain  row- 
charge  cratering  experience  in  a  weak, 
wet  medium.  In  addition,  this  experiment 
tested  techniques  for  connecting  a  row- 
charge  crater  to  an  existing  crater  and 
for  over -excavating  to  accept  throw out 
from  a  follow-on  connecting  row -charge 
crater. 


Project  Pre -Gondola  II  was  detonated 
at  Valley  County,  near  the  edge  of  the 
Fort  Peck  Reservoir  approximately  18  mi 
south  of  Glasgow,  Montana,  at 
0800  hr  (MDT),  28  June  1967  (see  Fig.  1). 
Coordinates  of  the  center  charge  were 
1V106°38»31”,  N47°55'51’'.  The  orientation 
of  the  row  was  along  an  alignment  10® 

East  of  North  to  10®  West  of  South.  This 
alignment  extended  through  the  center  of 
the  Charlie  crater  (an  existing  crater 
created  during  the  Pre -Gondola  I  series). 

The  average  lip  crest-to-lip  crest 
dimensions  after  connection  to  the  Charlie 
crater  were  640  ft  X  280  ft.  Individual 
charge  yields,  depths,  spacing,  and  re¬ 
sulting  apparent  crater  dimensions  are 
given  in  Table  I. 


Table  I,  Pre-Gondola  II  charge  yields,  depths,  spacing,  and  resulting  apparent  crater 
dimensions. 


Charge 

Tons  NM^ 

Emplacement 

configuration 

Apparent 

crater  dimensions 
(ft) 

Width  W  Depth  D 

a  a 

Depth  of 
burst 

Spacing 

Charlie  (existing  crater)^ 

(19.62) 

(42.5) 

r» 

(160.8) 

(32.6) 

E  (connecting  charge) 

38.61 

59.7 

105.5^ 

206.5 

55.5 

F 

19.70 

49.4 

7  9.8 

152.5 

37.5 

G  (center  charge)* 

19.55 

48.8 

79.9 

164.0 

36.9 

H  (over-excavating  charge) 

39.56 

59.9 

7  9.9 

r7  n  r\ 

214.5 

57.0 

I 

20.00 

48.8 

7  9.  9 

173.0 

33.5 

^Spherical  charges  of  liquid  explosive  nitromethane  (CH3NO2).  Actual  total  charge 
weight  was  137.42  tons. 

^Dimensions  for  Charlie  crater  are  those  existing  prior  to  row -charge  event. 


^Distance  between  this  charge  and  charge  on  line  above. 
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Fig.  1.  Map  showing  environs  of  Pre-Gondola  II  and  associated' 
events. 
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DESCRIPTION  OF  EVENTS  ASSOCIATED 
WITH  PRE -GONDOLA  II 

Prior  to  the  row -charge  detonation, 
two  1000 -lb  shots  were  fired  to  investigate 
the  effects  of  small -magnitude  decoupling 
factors.  These  events  were  SD-1  (fully 

Charge  Cavity  volume 

Event  (tons  NM)  explosive  volume 

SD-1  0.50  1.0 

SD-2  0.50  2.0 

DESCRIPTION  OF  OTHER  EVENTS  RE¬ 
LATED  TO  PRE -GONDOLA  II 

In  addition  to  the  above  events,  con¬ 
siderable  additional  analysis  has  been 
made  of  seismic  records  taken  during 
Pre -Gondola  I,  located  at  the  same  site 
as  the  Pre -Gondola  II  events.  This 


analysis 

is  reported  here. 

The  work  was 

Charge 

DOB 

Event 

(tons  NM) 

(ft) 

Charlie 

19.62 

42.49 

Bravo 

19.36 

46.25 

Alfa 

20.35 

52.71 

Delta 

20.24 

56.87 

OBJECTIVES 

In  order  to  provide  a  knowledgeable 
basis  to  plan  for  possible  large  scale  geo- 
nuclear  projects  in  the  future,  it  is  nec¬ 
essary  to  gather  information  on  the 
characteristic  behavior  of  various  earth 
materials  at  the  shot  point  and  in  the 
vicinity  of  subsurface  explosions.  The 
Pre -Gondola  series  of  explosions  offers 
an  opportunity  to  study  the  effect  of  several 
important  parameters  on  the  transmitted 


coupled)  and  SD-2  (decoupled, 

V  =2.0).  These  events  were  both 
exp 

detonated  on  14  Jme  1968  at  the  same  site 
as  the  140 -ton  event.  Cratering  data  for 
these  two  events  are  summarized  as 
follows: 


DOB 

Apparent  crater  dimensions 
(ft) 

(ft) 

Width  Wa 

Depth  Da 

17.3 

50.2 

10.5 

17.7 

47.4 

10.5 

originally  reported  in  PNE-1105.^  Since 
the  writing  of  that  report  additional 
analysis  techniques  have  become  available, 
and  it  is  imperative  that  for  a  complete 
\inder standing  of  the  results  of  Pre- 
Gondola  II,  these  data  and  results  be  in¬ 
cluded  at  this  time.  A  summary  of  the 
Pre -Gondola  I  cratering  data  is  as  follows: 


Apparent  crater  dimensions 
_ (ft) _ 


Width  W^ 

cx 

Depth 

160.8 

32.6 

157.0 

29.5 

152.2 

32.1 

130.2 

25.2 

seismic  signals. 

The  variety 

of  yields,  depths  of  burst,  and  shot  config¬ 
urations,  all  being  detonated  at 
essentially  the  same  geographical  location, 
makes  it  possible  to  make  repeated  mea¬ 
surements  at  the  same  locations.  This 
process  means  that  the  effect  of  trans¬ 
mission  path  can,  to  a  large  extent,  be 
normalized  out  of  the  data  at  any  given 
station.  (This  is  not  entirely  true  since 
the  ’’transmissibility'*  of  the  transmission 
path  is  frequency-dependent  and,  therefore. 
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may  "change"  if  the  source  fvmction  is 
altered  drastically  with  respect  to 
frequency — see  Appendix  A.) 

The  explosions  were  detonated  in  a  wet 
clay-shale  media  which  is  characteristic 
of  several  locations  of  interest  to  the  nu¬ 
clear  cratering  program.  This  region  of 
northeastern  Montana  is  also  of  interest 
because  it  appears  to  be  relatively  simple 
from  a  physical  properties  standpoint. 
Flat-lying,  extensive,  thick  layers  of 
sedimentary  materials  exist  to  depths  of 
over  8,000  ft  and  distances  of  over  100  km. 

In  view  of  the  above,  the  following  objectives 
were  chosen  for  this  study: 

1.  Dependence  of  velocity  time  history 
of  ground  motion  upon  yield,  depth  of 
burst,  and  various  shot  configurations 
(single  charge,  row -charge,  and  small 
decoupling  ratios). 

2.  Dependence  of  peak  velocity  spectra 
and  power  spectral  densities  upon  yield, 
depth  of  burst,  and  shot  configuration. 

3.  Identification  of  those  properties  of 
explosively  induced  seismic  signals  which 
scale  uniformly  and  consistently  over  large 
ranges  of  yield,  depth  of  burst,  and  varia¬ 
tions  in  shot  configurations. 

4.  Empirical  measurements  of  seismic 
motion  from  cratering  explosions  in  an  ex¬ 
tensive  water -satxirated  clay-shale  mate¬ 
rial. 

5.  A  study  of  the  attenuation  with  dis¬ 
tance  of  peak  signal  amplitudes,  velocity 
spectra,  and  power  spectral  density  in  a 
simple  geology. 

BACKGROUND 

A  word  of  explanation  is  necessary  at 
this  point.  Although  this  report  is  tech¬ 
nically  a  report  on  Pre -Gondola  II,  it  also 


contains  much  information  gathered  during 
Pre -Gondola  I.  The  reasons  for  this  are 
two-fold:  (1)  the  events  from  the  separate 
programs  are  meaningful  to  the  fullest 
extent  only  when  correlated  and  presented 
as  a  unit,  and  (2)  the  power  spectral 
density  analysis  technique  (PSD)  was  not 
available  when  the  Pre -Gondola  I  work  was 
reported.  It  was  the  conclusion  of  that 
report  (PNE-1105)  that,  when  the  analysis 
technique  did  become  available,  the  data 
should  be  reported.  It  is  included  in  this 
report  to  simplify  the  task  of  both  the 
reader  and  writer. 

Figure  1  shows  the  location  of  the  Pre- 
Gondola  experiments  in  northeastern 
Montana.  The  ideal  geology  for  seismic 
propagation  studies  would  be  a  simple  two- 
layer  medium,  flat,  and  of  infinite  extent 
in  the  horizontal  directions.  This  is  the 
simplest  geometry  that  would  allow  for 
propagation  of  surface  waves,  direct  p  and 
s  arrivals  and  refracted  arrivals.  The 
author  is  unaware  of  the  existence  any¬ 
where  of  such  a  geologic  setting  of 
appropriate  dimensions. 

The  geology  of  northeastern  Montana  is 
however  relatively  simple.  Massive 
sedimentary  beds  make  up  the  surface 
geology.  These  vary  in  thickness  up  to 
2000  ft  and  comprise  in  total  thickness 
about  5000-8000  ft.  Structural  features 
such  as  the  Willi ston  Basin  are  of  the  order 
of  himdreds  of  miles  across  and  of  relative¬ 
ly  shallow  relief,  resulting  in  relatively 
flat  bedding  planes. 

Previous  work  on  seismic  attenuation 
has  been  primarily  concerned  with  peak 
amplitudes  of  motion  in  the  wave  trains 
with  little  or  no  regard  for  the  total  fre¬ 
quency  content  of  the  motion. The  rea¬ 
sons  for  this  condition  have  been  twofold. 
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Until  recently  there  has  been  a  lack  of 
good  instrumentation  and  analysis  capa¬ 
bilities  to  do  this  work.  More  important, 
most  of  the  experience  to  date  has  been  in 
such  complicated  geologic  settings  that 
meaningful  frequency  analysis  has  been 
almost  impossible.  The  relatively  simple 
geology  offers  an  opportunity  to  investigate 
both  the  frequency  content  and  the  amplitude 
as  a  function  of  distance  from  the  source 
and  for  variations  in  yield  and  depth  of 
burst  of  the  source. 

Most  explosions  with  yields  greater 
than  5  tons  have  fallen  into  three  categories: 

1.  Industrial  explosions  of  all  kinds 
with  usually  little  or  no  instrumenta¬ 
tion. 

2.  Well -instrumented,  completely  con¬ 
tained  single -event  explosions. 

3.  Well-instrumented  cratering  exper¬ 
iments. 

For  the  most  part  the  larger  cratering 
experiments  prior  to  Pre -Gondola  have 
been  single  events.  Where  there  have 
been  several  events,  the  layout  and 
seismic  instrumentation  offered  no 
opportionity  to  study  coupling  of  energy 
into  the  seismic  signal. 

A  requirement  for  a  dynamic  cratering 
condition  (as  opposed  to  a  subsidence 
crater  )  is  that  the  cavity  communicates 
with  the  atmosphere  before  the  cavity 
growth  has  stopped.  This  means  that  the 


overburden  must  fail  in  tension,  and  the 
thus -isolated  blocks  of  overburden  will  be 
ejected  along  a  ballistic  trajectory. 
Seismically,  this  differs  from  the  com¬ 
pletely  contained  situation  by  allowing  the 
cavity  pressure  to  drop  at  early  times, 
thus  preventing  the  coupling  of  some 
fraction  of  the  energy  into  the  seismic 
signal  at  late  times.  The  possible  de¬ 
crease  in  seismic  energy  is  an  unknown, 
and  previous  estimates  of  seismic  motions 
for  cratering  shots  have  been  based  on  re¬ 
sults  for  the  completely  contained  case. 

If  there  is  an  effect  due  to  cratering  it 
could  be  proportional  to  a  function  of 
either  the  actual  depth  or  the  scaled  depth 
(D/W^^^*"^).  In  either  case  some  effect 
would  be  noticed  for  a  shallow,  low-yield 
explosion.  If  coupling  is  a  function  of 
actual  depth  (in  other  words,  a  function  of 
the  total  time  before  breakout),  one  would 
expect  to  see  a  less  marked  effect  for 
deeper,  high-yield  events.  If  the  effect 
is  a  function  of  the  scaled  depth,  then  it 
should  be  the  same  for  all  yields.  Pre- 
Gondola  II  and  associated  events  offer  a 
suitable  variety  of  events  to  study  these 
possibilities. 


A  subsidence  crater  occurs  when  an 
explosively  created  subsurface  cavity 
collapses  resulting  in  a  surface  depres¬ 
sion  without  any  direct  cavity -atmosphere 
interaction. 


Procedure 


STATION  SELECTION 

The  number  and  location  of  seismic 
stations  for  Pre -Gondola  II  were  essentially 
the  same  as  for  Pre -Gondola  I,  since  the 
criteria  were  identical.^  Figure  2  is  a 
map  showing  the  location  of  stations  for 
the  events  reported  herein.  All  stations, 
their  locations  with  respect  to  SGZ,  and 
the  events  during  which  they  were  occupied 
are  given  in  Table  II  along  with  other  infor¬ 
mation.  This  coverage  is  sufficient  to 
provide  the  following: 


1.  Five  stations  on  a  radius  for  two 
different  yields  (one  of  them  a  row 
charge) 

2.  One  station  repeated  for  a  single 
yield  at  four  depths  of  burst 

3.  Two  stations  repeated  for  the  de¬ 
coupling  sequence 

4.  One  station  repeated  for  three  yields 
at  optimum  depth  of  burst 

5.  Stations  90  deg  apart  with  respect 
to  SGZ  to  study  the  asymmetry  of 
seismic  radiation  from  a  row 
charge. 


Table  II.  Events  reported  and  data  analyses  used  for  velocity  geophone  measurements. 


Event 

Station 

Station  location 
Distance  Azimuth 

Quality 
of  data 

Analysis  tech¬ 
niques  used^ 

Pre -Gondola  I 

Alfa 

2BN 

11.65 

80E  of  N 

Excellent 

VA,  PSD 

Charlie 

2BN 

11.65 

8°E  of  N 

Excellent 

VA,  PSD 

Delta 

2BN 

11.65 

8‘>E  of  N 

Excellent 

VA,  PSD 

Bravo 

2BN 

11.65 

S^E  of  N 

Excellent 

VA,  PSD 

Bravo 

3N 

29.3 

9°E  of  N 

Poor 

VA,  PSD 

Bravo 

4N 

52.0 

8“E  of  N 

Good 

VA,  PSD 

Bravo 

5N 

95.9 

6®E  of  N 

Excellent 

VA,  PSD 

Pre -Gondola  II 

SD-l 

2BN 

11.65 

8»E  of  N 

Fair 

VA,  CFPV,  PSD 

SD-1 

6W 

6.39 

90®W  of  N 

Good 

VA,  CFPV,  PSD 

SD-2 

2BN 

11.65 

8°E  of  N 

Fair 

VA,  CFPV,  PSD 

SD-2 

6W 

6.39 

90°W  of  N 

Good 

VA,  CFPV,  PSD 

Row  charge 

2BN 

11.65 

8'>E  of  N 

Excellent 

VA,  CFPV,  PSD 

(E.  F,  G,  H,  I) 

3N 

29.3 

9°E  of  N 

Excellent 

VA,  CFPV,  PSD 

4N 

52.0 

8“E  of  N 

Excellent 

VA,  CFPV,  PSD 

5N 

95.9 

6°E  of  N 

Excellent 

VA,  CFPV,  PSD 

IN 

'-5.0 

8‘>E  of  N 

Very  poor 

None 

6W 

6.39 

90°W  of  N 

Very  poor 

None 

^VA  =  vector  addition,  PSD  =  power  spectral  density,  CFPV  =  comb  filter  peak 
velocity. 
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Fig.  2.  Map  showing  location  of  free-field  seismic  stations  for  Pre-Gondola 
II  140 -ton  row  charge.  (Station  2BN  was  operated  for  Alfa,  Charlie, 
Delta,  and  Bravo  Events,  and  Stations  3N,  4N,  and  5N  were  occupied 
on  Bravo  Event  d\xring  Pre -Gondola  series;  Stations  2BN  and  6W  were 
occupied  on  decoupling  series  during  Pre -Gondola  IL) 
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SEISMIC  INSTRUMENTATION 

Ground  motions  were  measured  with 
modified  HS-10-1  velocity  geophones.  The 
HS-10-1  is  a  moving  coil  instrument  with 
a  l“Sec  natural  period;  with  70%  damping 
it  has  a  flat  response  curve  above  1  Hz 
(cycles  per  second).  The  sensitivity  is 
determined  by  the  coil  impedance  (Fig.  3). 
To  improve  performance  and  reliability 
and  to  reduce  instrument  drift,  adjustments 
were  made  in  the  field  and  all  instruments 
were  monitored  after  emplacement. 
Emplacement  was  accomplished  about 
1  wk  prior  to  the  event  date,  and  the  final 
monitor  check  was  made  on  the  evening 
prior  to  the  event  date.  The  critical 


parameters  are  sensitivity,  natural  period, 
and  damping  ratios. 

Each  station  consisted  of  one  vertical 
and  two  horizontal  instruments,  the  latter 
oriented  radially  and  transversely  to  the 
shot  point.  This  comprises  the  standard 
rectangular  system,  giving  motion  in  the 
three  mutually  orthogonal  directions 
necessary  to  define  translational  motion. 

As  a  matter  of  definition,  for  the  rest  of 
this  report  the  following  designations  will 
be  used:  x  refers  to  radial  motion,  y 
refers  to  transverse  motion,  z  refers  to 
vertical  motion,  and  Res  is  the  vector  sum 
(n/x^  +  y^+  z^)  for  velocity  and  an 
arithmetic  sum  (x  +  y  +  z)  for  energy. 

Isolation  from  wind  noise  and  good 
mechanical  coupling  with  the  earth  are 


Frequency  —  Hz 


Fig.  3.  Response  of  HS-lO-l  velocity  geophone  at  35  and  70%  damping. 
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accomplished  by  burying  the  instruments. 
To  meet  the  requirements  of  instrument 
protection,  ease  of  operation  and  emplace¬ 
ment,  quality  of  orientation,  etc.,  an 
instrument  canister  was  developed  to 
maintain  the  instruments  in  the  mutually 
orthogonal  relationship.  The  water-tight 
canister  is  then  easily  buried  in  a  12 -in. 
diameter,  40 -in.  hole  and  oriented  as  a 
unit.  Tamped  sand  was  used  as  a  backfill. 
Figure  4  illustrates  the  details  of  the 
instrument  package. 

RECORDING  AND  TIMING 

Magnetic  tape  recorders  utilizing  an 
FM  electronics  system,  a  four-track  head. 


and  a  1/4 -in.  loop -tape  cartridge  with  a 
5 -1/2 -min  capacity  were  used  on  all 
measurements.  One  of  these  recorders 
is  shown  in  Fig.  5  and  partly  dismantled 
in  Fig.  6.  The  recorder  shown  was  used 
on  the  Pre -Gondola  I  Events  Alfa,  Bravo, 
Charlie,  and  Delta.  For  the  Pre-Gondola 
II  Events  certain  modifications  were  made 
on  the  power  supplies  and  VCO’s  to  im¬ 
prove  their  temperature  and  time 
stability.  These  modifications  did  not 
significantly  alter  the  appearance  of  the 
recorders. 

Each  geophone  (component)  was  recorded 
on  one  track  of  the  magnetic  tape  and  an 
oscillator -controlled  timing  signal  was 
recorded  on  the  fourth  track.  An 


Fig.  4.  Three -component  instrument  package  during  set-up  and  prior  to  in¬ 
sertion  into  outer  canister  (canister  is  partly  visible  behind  metal 
drum  on  right). 
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Fig.  5.  Top  view  of  tape  recorder  used  for  recording  geophone  signals  (tape  cartridge 
is  shown  in  place  at  right). 


approximate  zero  time  signal  (square 
pulse)  was  manually  superimposed  on  the 
timing  channel  by  the  station  operator  who 
was  in  radio  contact  with  the  firing  control 
center. 

The  cross  timing  accuracy  between 
stations  is  determined  by  the  variation  in 
station  operator  reaction  time.  This 
variation  could  cause  an  uncertainty  of  as 
much  as  1  sec  in  actual  zero  time.  This 
uncertainty  is  considered  too  large  to 
allow  meaningful  analysis  of  the  transit 
times  of  various  wave  phases,  etc. 

INFORMATION  RETRIEVAL  AND 
REDUCTION 

Instrument  responses  to  step  functions 
were  recorded  prior  to  events.  Recorder 


calibration  signals  (sine  waves  and  ramps) 
were  recorded  on  the  "event  tapes"  before 
and  after  the  event.  All  recording  equip¬ 
ment  was  turned  on  and  running  at  minus 
30  sec.  A  positive  dc  voltage  of  about 
1-sec  duration  was  manually  superim¬ 
posed  at  zero  time  on  the  timing  channel. 

Subsequent  to  the  event,  magnetic  tape 
records  from  all  stations  were  redubbed 
to  a  1  -in.  tape  to  make  them  compatible 
with  the  digital  analysis  system.  Records 
were  digitized  with  an  Astro-Data  analog- 
to-digital  translator.  Analysis  of 
digitized  records  was  completed  on 
IBM  7094  and  CDC  3600  computers.  The 
digitized  record  contained  400  data  points 
per  sec  for  each  component. 


DATA  ANALYSIS 

The  most  used  method  of  analyzing 
seismic  motion  data  of  the  type  with 
which  we  are  concerned  has,  in  the  past, 
been  to  measure  the  peak  amplitude 
attained  regardless  of  where  in  the  record 
it  occurred.  Such  an  analysis  is  crude  to 
be  sure,  but  little  else  could  be  done  due 
to  the  high  cost  of  hand  analysis  on  the 
many  records  that  had  to  be  studied. 
Digital  computers  have  become  common 
enough  to  reduce  the  total  cost  and  more 


sophisticated  analysis  techniques  are  now 
possible. 

It  is  of  interest  at  this  time  to  compare 
this  method  to  newer  ones  and  to  obtain 
some  measure  of  the  validity  of  the  old 
method.  Accordingly  peak  motions  were 
measured  from  the  records.  Also  mea¬ 
sured  were  the  time  of  arrival  of  that 
motion  after  the  first  motion  and  the 
approximate  frequency  of  the  motion. 

The  latter  is  calculated  by  assuming  the 
half  wave  to  be  of  harmonic  nature;  i.  e.. 


f  =  (1/2)  (1/At)  where  At  is  the  time  be¬ 
tween  zero  crossings  just  prior  to  and 
following  the  peak  motion. 

One  of  the  newer  techniques  for 
analyzing  a  time  varying  signal  is  to  filter 
it  into  its  spectral  components  to  mea¬ 
sure  the  peak  amplitudes  in  each  compo¬ 
nent  generating  a  velocity  spectra.  This 
method  is  called  comb  filter  peak 
velocity  (CFPV)  and  is  described  in 
Appendix  B. 

The  actual  equipment  used  in  this 
analysis  was  band  pass  filters.  The 
signal  ratios  at  1  octave  from  the  center 
frequency  were  about  15  dB  for  power 
(voltage^^/ voltage^^^  «  6).  Center  fre¬ 
quencies  were  set  in  a  semi -systematic 
pattern  to  cover  the  range  of  interest. 

They  were  0.75,  1.0,  1.5,  2.0,  3.0,  4.0, 
5.0,  6.0,  8.0,  10.0,  12.0,  and  15.0  Hz, 
respectively. 

Comb  filtering  was  planned  for  the 
records  obtained  in  the  140-ton  row 
charge,  as  well  as  the  SD-1  and  SD-2 
Events.  The  peak  amplitudes  for  each 
frequency  from  a  given  geophone  were 
plotted  against  frequency  to  illustrate  the 
spectral  content  of  the  signal.  Comparison 
of  the  plots  of  the  three  space  coordinates, 
the  different  stations,  and  the  several 
events  provides  an  understanding  of  how 
the  spectral  content  varies  with  distance, 
space  direction,  source  configuration  (for 
the  SD-1  Event  compared  to  the  SD-2)  and 
for  yield  (the  SD-1  Event  compared  to  the 
140 -ton  row  charge).  Attenuation  for 
various  frequencies  can  be  determined 
from  a  fit  to  the  data.  The  equation  is 


A^  =  Amplitude  (f)  = 

where  a  is  the  attenuation  slope. 

A  comparison  of  these  results  to  those 
obtained  by  the  power  spectral  density 
method  of  analysis  (PSD)  provides  a  basis 
for  choosing  the  most  consistent  and 
appropriate  method  of  the  two  for  analyz- 
ing  ground  motions  from  under groxand 
explosions. 

The  three  rectangular  components  are 
sufficient  to  describe  the  translational 
motion  of  a  point.  However  to  gain  a 
knowledge  of  the  absolute  amplitude  of 
motion  requires  a  transformation  to 
spherical  coordinates.  The  theory  and 
computer  code  for  the  transformation  is 
described  in  Appendix  B. 

The  code  calculates  x,  y,  z,  (Res)  , 
|Res|,  6,  and  c|)  as  functions  of  time.  It 
is  called  (surprisingly)  VECTOR. 

Another  method  of  analyzing  a  signal 
to  obtain  a  measure  of  its  spectral  con¬ 
tent  and  total  energy  is  to  determine  the 
PSD  and  the  integral  of  (PSD)df  or  E. 

This  method  is  described  in  Appendix  C 
and  Ref.  5.  The  quantities  produced  by 
the  PSD  code  for  a  given  time  history  of 
velocity  are  PSD(f),  or  density  of  power 
as  a  function  of  frequency  for  three 
components,  a  resultant  E,  or  total 
energy  in  a  component,  and  J^E  or  total 
station  energy. 

The  PSD*s  for  all  signals  in  this  report 
were  to  be  calculated  and  comparisons 
made  to  determine  the  variance  with  yield, 
depth  of  burst,  source  configuration,  and 
radius  and  to  make  selected  comparisons 
with  the  CFPV’s  for  some  events  and 
stations. 


Results 


DATA  RECOVERY  AND  REDUCTION 

Table  II  lists  the  events  included  in 
this  report  along  with  the  various  seisnaic 
stations,  their  locations  with  respect  to 
SGZ,  the  quality  of  the  data  recorded, 
and  the  analysis  techniques  used  on  each 
station  recording. 

The  failure  of  the  primary  recorder 
is  the  reason  for  the  poor  quality  on 
Station  3N  for  Bravo,  A  backup  recorder 
set  at  low  gain  recorded  the  event  but  the 
signal-to -electronic  noise  ratio  is  small. 
On  the  140-ton  row  charge  the  motions 
at  Stations  6W  and  IN  exceeded  the 
mechanical  limits  of  the  geophones, 
making  the  data  totally  unreliable.  No 
results  for  these  two  stations  on  this 
event  will  be  included.  The  small  motions 
from  Events  SD-1  and  SD-2  resulted  in  a 
small  signal-to-background  seismic 
noise,  reducing  the  quality  of  these  signals 
somewhat. 

The  method  of  data  reduction  requires 
redubbing  which  introduces  additional 
noise,  usually  random  or  white,  but 
sometimes  peculiar  to  a  given  machine 
as  a  result  of  mechanical  flutter  in  the 
tape  drive,  etc.  This  is  apparently  the 
cause  of  spikes  in  the  power  spectral 
density  curves  at  9-1/2  and  15  Hz  in  all  of 
the  SD-1,  SD-2,  and  140-ton  row-charge 
Events  from  Pre -Gondola  II.  The  large 
spike  in  the  PSD  for  all  components  of 
Station  3N  on  Bravo  is  not  present  for  the 
140 -ton  Event  at  3N  nor  is  it  present  on 
any  other  station.  It  is  assumed  that  this 
is  due  to  the  tape  drive  of  the  event  re¬ 
corder  especially  since  this  was  the  low 
gain  recorder.  There  is,  however,  some 


doubt  about  the  validity  of  this  assumption 
and  both  possibilities  are  reported  in  the 
peak  amplitude  values.  It  has  little  effect 
upon  the  total  energy  E  (see  Appendix  C) 
in  the  signal. 

PEAK  VELOCITY  AND  VECTOR 
ADDITION 

Table  III  lists  the  peak  velocity  of  the 
components  and  the  peak  vector  resultant, 
the  time  of  arrival  of  that  peak  after  first 
motion,  and  the  apparent  frequency  of 
that  motion  for  all  stations  and  all  events. 
Peak  velocities  for  the  stations  recording 
for  Bravo  and  the  140 -ton  row  charge  are 
plotted  in  Figs.  7  and  8  as  a  function  of 
distance  from  SGZ. 

COMB  FILTER  PEAK  VELOCITY 

Only  the  records  from  SD-1,  SD-2, 
and  the  140 -ton  row  charge  were  analyzed 
in  this  manner.  Table  IV  contains  these 
data.  A  better  concept  of  these  data  is 
gained  when  they  are  put  in  graphical 
form.  Figures  9  and  10  show  the  results 
for  Station  2BN  on  SD-1  and  SD-2,  re¬ 
spectively,  and  Figs.  11,  12,  and  13  com¬ 
pare  the  three  separate  space  coordinates 
for  the  two  events  to  illustrate  the  differ¬ 
ences  due  to  the  small  decoupling  factor. 
Figures  14  through  18  show  results  for 
Station  6W  for  the  same  events  in  the 
same  sequence. 

The  CFPV -frequency  curves  for  each 
station  for  the  140 -ton  row -charge  are 
shown  in  Figs.  19  through  22,  while  the 
individual  components  for  all  stations  are 
compared  in  Figs.  23  through  25.  The 
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Table  III.  Velocity  data — all  stations. 


Event 

Station 

Component 

Peak  velocity 
(mm/ sec) 

Time  of  arrival 
after  first  motion 
(sec) 

Frequency 
of  motion 
(Hz) 

SD-1 

2BN 

X 

0.0385 

11.2 

1.75 

SD-1 

2BN 

y 

-'0.027 

11.6 

1.75 

SD-1 

2BN 

z 

0.0338 

0.8 

6.8 

SD-1 

2BN 

Res 

0.042 

11.6 

- 

SD-2 

6W 

X 

0.108 

7.6 

2.4 

SD-2 

6W 

y 

0.096 

9.4 

3.7 

SD-2 

6W 

z 

0.096 

9.97 

8.0 

SD-2 

6W 

Res 

0.113 

7.6 

- 

SD-2 

2BN 

X 

0.0545 

11.9 

1.75 

SD-2 

2BN 

y 

-0.045 

12.5 

1.75 

SD-2 

2BN 

z 

0.0432 

1.4 

7.6 

SD-2 

2BN 

Res 

0.0687 

12.6 

- 

SD-1 

6W 

X 

0.055 

7.5 

2.5 

SD-1 

6W 

y 

0.072 

7.5 

3.4 

SD-1 

6W 

z 

0.099 

1.0 

8.8 

SD-1 

6W 

Res 

0.099 

1.0 

- 

Bravo 

2BN 

X 

0.76 

12.4 

1.75 

Bravo 

2BN 

y 

0.34 

12.4 

1.2 

Bravo 

2BN 

z 

1.20 

2.4 

4.5 

Bravo 

2BN 

Res 

1.29 

2.4 

- 

Bravo 

3N 

X 

0.16 

? 

1.5 

Bravo 

3N 

y 

0.10 

? 

1.5 

Bravo 

3N 

z 

0.125 

? 

2.5 

Bravo 

3N 

Res 

0.220 

? 

- 

Bravo 

4N 

X 

0.0207 

12.65 

3.0 

Bravo 

4N 

y 

0.0142 

14.7 

3.0 

Bravo 

4N 

z 

0.0257 

13.5 

2.5 

Bravo 

4N 

Res 

0.0306 

13.5 

- 

Bravo 

5N 

X 

0.024  ±  10% 

5.48 

7.0 

Bravo 

5N 

y 

0.014  ±  10% 

5.53 

8.0 

Bravo 

5N 

z 

0.08  ±  10% 

4.2 

8.0 

Bravo 

5N 

Res 

0.08  ±  10% 

4.2 

- 

Alfa 

2BN 

X 

1.10 

2.32 

6.0 

Alfa 

2BN 

y 

0.45 

11.9 

2.0 

Alfa 

2BN 

z 

1.94 

2.25 

5.0 

Alfa 

2BN 

Res 

2.09 

2.25 

- 

Charlie 

2BN 

X 

0.75 

3.2,  11.2 

5.0,  1.5 

Charlie 

2BN 

y 

0.28 

12.8 

1.5 

Charlie 

2BN 

z 

1.38 

3.1 

5.0 

Charlie 

2BN 

Res 

1.46 

3.1 

Delta 

2BN 

X 

1.72 

11.8 

1.75 

Delta 

2BN 

y 

0.72 

13.5 

1.75 

Delta 

2BN 

z 

2.20 

2.26 

6.00 

Delta 

2BN 

Res 

2.45 

2.34 

- 

Row  charge 

2BN 

X 

4.4 

10.7 

2.0 

Row  charge 

2BN 

y 

2.8 

11.7 

0.8 

Row  charge 

2BN 

z 

5.8 

2.21 

4.8 

Row  charge 

2BN 

Res 

6.3 

2.21 

- 

Row  charge 

3N 

X 

0.251 

10.5 

2.0 

Row  charge 

3N 

y 

0.209 

,  20.8 

1.75 

Row  charge 

3N 

z 

0.655 

6.8  ' 

3.0 

Row  charge 

3N 

Res 

0.663 

6.8 

- 

Row  charge 

4N 

X 

0.193 

12.6 

3.1 

Row  charge 

4N 

y 

0.14 

14.2 

3.0 

Row  charge 

4N 

z 

0.32 

12.0 

2.7 

Row  charge 

4N 

Res 

0.344 

12.2 

Row  charge 

5N 

X 

0.097 

5.4 

4.0 

Row  charge 

5N 

y 

0.046 

12.9 

3.3 

Row  charge 

5N 

z 

0.197 

4.0 

7.0 

Row  charge 

5N 

Res 

0.197 

4.0 

- 

-14- 


mm /sec 


100 


Radius  —  km 


Fig.  7.  Peak  velocities  for  Pre -Gondola  I 
Bravo  Event,  as  a  function  of  dis¬ 
tance  (large  value  at  96  km  is  due 
to  single  pulse  of  very  high 
amplitude  in  velocity -time  signa¬ 
ture;  rest  of  signature  is  of  rela¬ 
tively  low  amplitude). 


Radius  —  km 

Fig.  8.  Peak  velocities  for  Pre -Gondola  II 
140 -ton  row  charge  as  a  function 
of  distance. 


Table  IV.  Comb  filter  peak  velocity  data  for  Pre -Gondola  II. 


Event 

Station 

Component 

0.75  Hz 

l.OHz 

1.5  Hz 

2.0  Hz 

3.0  Hz 

4.0  Hz 

5.0  Hz 

6.0  Hz 

8.0  Hz 

10.0  Hz 

12.0  Hz 

15,0  Hz 

140 -ton 

2BN 

X 

1.72 

2.87 

2.87 

2.10 

0.86 

1.62 

2.10 

1.91 

1.53 

0.96 

0.48 

0.19 

row  charge 

y 

0.67 

0.78 

1.34 

1.23 

0.56 

0.33 

0.67 

0.56 

0.47 

0.33 

0.23 

0.12 

z 

0.97 

1.84 

2.16 

1.30 

4.00 

5.40 

5.40 

4.64 

2.38 

1.73 

0.97 

0.54 

3N 

X 

0.23 

0,25 

0.18 

0.27 

0.21 

0.20 

0.14 

0.083 

0.062 

0.041 

0.021 

- 

y 

0.094 

0.19 

0.17 

0.15 

0.11 

0.12 

0.11 

0.11 

0.047 

0.024 

0.012 

■  - 

z 

0.19 

0.31 

0.20 

0.36 

0.63 

0.45 

0.22 

0.18 

0.13 

0.078 

0.045 

0.033 

4N 

X 

0.052 

0.042 

0.057 

0.12 

0.18 

0.15 

0.083 

0.052 

0.037 

0.026 

0.021 

0.010 

y 

0,029 

0.046 

0.032 

0.072 

0.13 

0.075 

0.046 

0.032 

0.023 

0,013 

0.009 

- 

z 

0.065 

0.054 

0.081 

0.20 

0.35 

0.23 

0.11 

0.065 

0.051 

0.070 

0.076 

0.085 

5N 

X 

0.022 

0.011 

0.022 

0.022 

0.039 

0.044 

0.039 

0.033 

0.033 

0.028 

0.022 

- 

y 

0.013 

0.007 

- 

0.013 

0.027 

0.013 

0.040 

0.034 

0.013 

- 

- 

z 

- 

- 

0.011 

0.045 

0.068 

0.057 

0.085 

0.11 

0.14 

0.13 

0.13 

0.068 

SD-1 

2BN 

X 

0.008 

0.010 

0.031 

0.036 

0.019 

0.013 

0.008 

0.010 

0.013 

0.008 

- 

_ 

y 

0.005 

0.005 

0.014 

0.014 

0,011 

0.007 

0.005 

0.011 

0.009 

0.005 

- 

- 

z 

0.005 

0.005 

0.012 

0.012 

0.010 

0.014 

0.017 

0.022 

0.026 

0.019 

0.014 

0.009 

6W 

X 

0.012 

0.018 

0.039 

0.058 

0.066 

0.066 

0.048 

0.051 

0.042 

0,030 

0.027 

0.012 

y 

0.007 

- 

0.014 

0.018 

0.053 

0.064 

0.043 

0.021 

0.028 

0.018 

0.014 

- 

z 

0.008 

0.011 

0.039 

0.044 

0.044 

0.048 

0,033 

0.048 

0.069 

0.066 

0.066 

0,061 

SD-2 

2BN 

X 

0.008 

0.019 

0,050 

0.053 

0.029 

0.015 

0,015 

0.015 

0.017 

0.009 

- 

- 

y 

0.005 

0.009 

0.023 

0.027 

0.014 

0.014 

0.011 

0.007 

0.018 

0.005 

- 

- 

z 

0.005 

0.012 

0.026 

0.026 

0.014 

0.017 

0.  024 

0.034 

0.036 

0.034 

0.014 

- 

6W 

X 

0.012 

0.030 

0.060 

0.072 

0.072 

0.072 

0.042 

0.036 

0.042 

0.033 

0.030 

0.015 

y 

0.007 

0.018 

0.036 

0.064 

0.078 

0.050 

0.028 

0.025 

0.025 

0.021 

0.014 

- 

z 

0.008 

0.025 

0.050 

0.056 

0.061 

0.044 

0.033 

0.050 

0.066 

0.077 

0.069 

0.061 
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ig,  9.  Comb  filter  peak  velocity  as  a  function  of  frequency  for  Station  2BN  Pre- 
Gondola  II  SD-1  Event. 


Fig.  10.  Comb  filter  peak  velocity  as  a  fvinction  of  frequency  for  Station  2BN, 
Pre -Gondola  II  SD-2  Event. 
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rig.  13.  Comparison  of  comb  filter  peak  velocities  of  vertical  motion  at  Station 
2BN  for  SD-1  and  SD-2  Events. 


Fig.  14.  Comb  filter  peak  velocity  as  a  function  of  frequency  for  Station  6W, 
Pre -Gondola  II  SD-1  Event. 
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Fig.  15.  Comb  filter  peak  velocity  as  a  function  of  frequency  for  Station  6W, 
Pre “Gondola  II  SD-2  Event. 


Fig.  16.  Comparison  of  comb  filter  peak  velocities  of  radial  motion  at  Station  6W^ 
for  SD-1  and  SD-2  Events. 
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Fig.  19.  Comb  filter  peak  velocity  as  a  fiinction  of  frequency  for  Station  2BN, 
Pre-Gondola  II  140-ton  row  charge. 


Fig.  20.  Comb  filter  peak  velocity  as  a  function  of  frequency  for  Station  3N, 
Pre-Gondola  II  140-ton  row  charge. 
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Fig.  23.  Comb  filter  peak  velocity  of  radial  components  for  all  stations,  Pre -Gondola  II 
140 -ton  row  charge. 


attenuation  a,  (see  Appendix  B)  with 
distance  for  the  peak  amplitudes  from  the 
CFPV’s  is  shown  in  Fig.  26. 

POWER  SPECTRAL  DENSITY 

Table  V  lists  the  peak  density 
the  frequency  at  which  it  occurs,  and  the 


total  energies  (E  and^E)for  each  com¬ 
ponent  on  all  stations  and  events. 

Comparisons  of  the  PSD’s  for  each 
component  (including  the  resultant)  at 
Station  2BN  for  three  yields  are  shown  in 
Figs.  27  through  30.  Comparisons  of  the 
PSD’s  for  each  component  and  resultant 
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Frequency  —  Hz 

Fig.  25.  Comb  filter  peak  velocity  of  vertical  components  for  all  stations,  Pre -Gondola  II 
140 -ton  row  charge,  as  a  function  of  frequency.  (Note  that  at  low  frequencies 
attenuation  is  about  as  expected;  at  high  frequencies  the  short  duration -high 
amplitude  pulses  in  signals  on  Station  5N  considerably  alter  relationship  of  V 
vs  distance.) 


(Pre -Gondola  II). 


Table  V.  Power  spectral  density:  peak  density  values  and  total  energy 


Event 

Station 

Component 

Peak  density 

Value  Frequency®- 

(mm/sec)^  sec  (Hz) 

Total  energy 
(E  and 
(mm/sec)2 

Pre -Gondola  II 

2  BN 

X 

0.805 

1.17 

0. 

,633  E 

~  140 -ton  row 

y 

0.164 

1.50 

0. 

.250  E 

charge 

z 

0,255 

1.33 

0. 

.570  E 

Res 

1.20 

1.17 

1, 

,453  2E 

3N 

X 

2.26 

X  10“3 

1.87 

5.65 

X  10’^ 

y 

2.72 

X  10“3 

1.67 

2.89 

X  10"3 

z 

7.74 

X  10~3 

2.67 

8.06 

X  i0~l 

Res 

10,5 

X  10^3 

2.67 

16.60 

X  10"® 

4N 

X 

1.24 

X  10-3 

2.67 

1.67 

X  10-3 

y 

3.43 

X  10-4 

2.83 

0.473  X  10-3 

z 

2.10 

X  10"3 

2.83 

2.76 

X  10-3 

Res 

3.55 

X  10“'=^ 

2.83 

4.903  X  10-3 

5N 

X 

1.54 

X  10'4 

2.83 

2.83 

X  10-4 

y 

0.274  X  10"4 

4.67 

1.36 

X  10-4 

z 

2.57 

X  10‘4 

2.67 

3.51 

X  10-4 

Res 

4.30 

X  10-4 

2.67 

7.70 

X  10-4 

SD-1 

2BN 

X 

4.68 

X  10'5 

1.83 

4.4 

X  10"® 

y 

2.51 

X  10-5 

1.63 

4.1 

X  10"® 

z 

1.62 

X  10-5 

1.50 

4.1 

X  10"® 

Res 

7.42 

X  10-5 

1.83 

12.6 

X  10-5 

SD-2 

2  BN 

X 

11.40 

X  10-5 

1.83 

8.8 

X  10“® 

y 

5.12 

X  lO'S 

1.67 

7.1 

X  10“® 

z 

4.06 

X  10“5 

1.50 

7.7 

X  10-5 

Res 

17.4 

X  10'5 

1.83 

23.6 

X  10“5 

6W 

X 

2.492  X  10“4 

1.67 

5.12 

X  10-4 

y 

3.28 

X  10-4 

2.17 

4.24 

X  10-4 

z 

2.35 

X  10-4 

3.00 

3.84 

X  10-4 

Res 

6.35 

X  10-4 

2.17 

13.2 

X  10-4 

SD-1 

6W 

X 

0.64 

X  10-4 

2.17 

1.31 

X  10-4 

y 

1.12 

X  10-4 

3.67 

1.95 

X  10-4 

z 

1.29 

X  10-4 

3.00 

2.67 

X  10-4 

Res 

2.45 

X  10-4 

3.16 

5.93 

X  10-4 

Pre -Gondola  I 

2  BN 

X 

2.17 

X  10“2 

1.17 

2.40 

X  10"2 

Bravo 

y 

0.66 

X  lO'*^ 

1.17 

0.70 

X  10-2 

z 

0.97 

X  10'2 

1.33 

2.05 

X  10-2 

Res 

3.61 

X  10-2 

1.33 

5.15 

X  10-2 

3N 

X  ^ 

f  9.09 

U.06 

X  10-4 

X  10-4 

7.50 

1.50 

29.03 

X  10-45 

(4.37 

X  10"4 

7.50 

12.91 

y  ' 

(2.17 

X  10‘4 

1.83 

X  10"^^ 

(3.94 

13.43 

X  10-4 

X  10-4 

7.50 

2.67 

14.69 

X  10-4^ 

Res  1 

17.40 

9.57 

X  10-4 

X  10"4 

7.50 

1.50 

56.63 

1 

o 

T*H 

X 

4N 

X 

14.77 

X  10-6 

3,50 

26.95 

X  10"® 

y 

2.92 

X  10"® 

2.83 

6.89 

X  10"® 

z 

13.58 

X  i0"6 

2.83 

23.34 

X  10“® 

Res  28.07 

X  10-6 

2.83 

57.18 

X  10"® 

5N 

X 

7.00 

X  10~6 

1.67 

16.14 

X  10-® 

y 

2.75 

X  10^ 

5.00 

8.96 

X  10-® 

z 

8.81 

X  10"® 

6.50 

39.88 

X  10-5 

Res  13.96 

X  10"^ 

2.67 

64.98 

X  10"® 

Alfa 

2BN 

X 

7.36 

X  10-2 

1.33 

8.58 

X  10-2 

y 

z 

0.736 

2.46 

X  io"5 

X  10'^ 

1.87 

1.50 

1.47 

6.91 

X  10"2 

X  10"2 

Res  10.26 

X  10-2 

1.33 

16.96 

X  10"2 

Charlie 

2BN 

X 

2.45 

X  10"2 

1.33 

2.67 

X  10-2 

y 

0,535 

X  10-2 

1.67 

0.65 

X  10-2 

z 

1.16 

X  10“2 

1.3  3 

3.02 

X  10"2 

Res 

3.96 

X  10"2 

1.33 

6.34 

X  lO"'^ 

Delta 

2BN 

X 

9.53 

X  10"2  , 

1.33 

10.46 

X  10"2 

y 

2.29 

X  10“2 

1.67 

2.16 

X  10-2 

z 

3.36 

X  10-2 

1.33 

7.68 

X  10"^ 

Res  13.60 

X  10-2 

1.33 

20.30 

X  10"-^ 

^Frequency  at  whibh  peak  density  occurs. 

^Spike  at  7-1/2  Hz  is  on  all  three  channels  and  unexplained  at  this  time. 
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Fig.  27.  Radial  components  at  Station  2BN  for  three  yields. 
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Resultant  power  spectral  density 
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at  Stations  6W  and  2BN  for  SD-1  and  SD-2 
Events  are  shown  in  Figs.  31  through  38. 

Comparison  of  PSD's  for  each  compo¬ 
nent  at  various  distances  are  shown  in 
Figs.  39  through  42  for  the  Bravo  Event, 
and  in  Figs.  43  through  46  for  the  140 -ton 
row  charge.  Figures  47  through  50  show 
comparisons  of  the  PSD*s  of  the  separate 


components  at  Station  2BN  for  the  20 -ton 
Pre -Gondola  Events  Alfa,-  Bravo,  Charlie, 
and  Delta  to  show  the  effect  of  depth  of 
burst. 

The  peak  PSD’s  for  Bravo  and  the  140- 
ton  row  charge  are  plotted  against  distance 
in  Figs.  51  and  52,  respectively.  The 
total  energies,  E  and  ^E,  are  shown  in 


Frequency  —  Hz 


Fig.  32.  Transverse  components  at  Station  6W  for  two  source  configurations. 


Figs.  53  and  54  for  the  same  events.  In 
a  manner  similar  to  that  done  for  CFPV*s, 
a  plot  of  the  distance  attenuation  param¬ 
eter  Q:(f)  of  the  PSD  components  for  each 
frequency  is  given  in  Fig.  55  for  the  Bravo 


Event  and  in  Ffg.  56  for  the  140 -ton  row 
charge.  A  comparison  of  the  a(f)  of  the 
resultant  PSD(f)  for  the  two  events  is  given 
in  Fig.  57. 

A  plot  of  a(f)/2  for  the  resultant  PSD 
is  compared  to  the  CFPV^s  a{f)  in  Fig.  58 


to  compare  the  results  of  the  two  spectral 
analysis  methods.  Since  there  is  no  re¬ 
sultant  component  of  CFPV,  an  approxi¬ 
mation  was  made  by  calculating  an  a{i) 
from  plots  of  peak  velocity  for  a  given  fre¬ 
quency  regardless  of  which  component  it 
occurred  in. 


Figure  59  shows  the  total  station  energy 
at  Station  2BN  for  the  Pre -Gondola  I 
Events  as  a  function  of  depth  of  burst. 
Figure  60  shows  the  total  energy  at  Station 
2BN  normalized  to  depth  of  burst  and 
plotted  vs  yield  for  three  yields  and  six 
events. 


Fig.  33.  Vertical  components  at  Station  6W  for  two  source  configurations 


--SD-1 


Frequency  —  Hz 


Fig.  34.  Resultant  power  spectral  densities  at  Station  6W  for  two  source  configurations, 


In  order  to  give  the  reader  a  feeling  for  shown  have  been  redrawn  for  reproduction 
the  actual  transducer  signals,  several  are  purposes  from  the  originals  and  are  not 
shown  in  Figs.  61  and  62.  (The  traces  suitable  for  analysis.) 
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Fig.  35.  Radial  components  at  Station  2BN  for  two  source  configurations. 
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Fig.  37.  Vertical  components  at  Station  2BN  for  two  source  configurations. 


Frequency  “  Hz 


Fig.  39.  Radial  components  for  Pre-Gondola  I  Bravo  Event  at  various  distances. 
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Fig.  40.  Transverse  components  for  Pre-Gondola  I  ] 
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Fig.  41.  Vertical  components  for  Pre -Gondola  I  Bravo  Event  at  various  distances. 
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Fig.  43.  Radial  components  for  Pre -Gondola  II  140 -ton  row  charge  at  various 
distances. 
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Power  spectral  density 


Fig.  47.  Radial  components  at  Station  2BN  for  four  20 -ton  events  at  different  depths 
of  burst. 
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Fig,  51.  Peak  power  spectral  densities 
plotted  against  radius  for  Pre- 
Gondola  I  Bravo  Event  (a  comes 
from  the  equation:  PSD„(R) 

=  KR-“).  P 
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Fig.  52.  Peak  power  spectral  densities 
plotted  against  radius  for  Pre- 
Gondola  II  140 -ton  row  charge 
{a  comes  from  the  equation: 
PSD  (R)  =  KR“«), 
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Fig.  53,  Energy  in  seismic  signal  for 

components  and  total  energy  as  a 
function  of  radius  for  Pre- 
Gondola  I  Bravo  Event  {a  comes 
from  the  equation: 


Fig.  54.  Energy  in  seismic  signal  for 

components  and  total  energy  as  a 
ftinction  of  radius  for  Pre- 
Gondola  II  140 -ton  row  charge 
{a  comes  from  the  equation: 
2E(R)  =  KR-“). 
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Fig.  55,  Slope  of  power  spectral  density  attenuation  as  a  function  of  frequency  for 

Event  (i.  e.,  slope  a  comes  from  the  equation: 
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Fig,  56.  Slope  of  power  spectral  density  attenuation  as  a  fionction  of  frequency  for 
PSD'(R°f??K^R-«?f))°'^  charge  (i.  e.,  slope  a  comes  from  the  equation: 
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Fig.  57.  Comparison  of  attenuation  for  resultant  power  spectral  densities  plotted 
function  of  frequency  for  Pre -Gondola  I  and  IT. 
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Fig.  58.  Comparison  of  slopes  of  attenuation  of  comb  filter  peak  velocities  to  one -half 
of  slopes  of  attenuation  of  power  spectral  density  resultant  for  Pre -Gondola  II 
140 -ton  row  charge  as  a  function  of  frequency. 
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Fig.  59.  Total  station  energy  (^E)  at 

Station  2BN  (Pre -Gondola  I)  show¬ 
ing  transmitted  energy  dependence 
upon  actual  depth  of  bTirst  (Bravo 
datum  is  neglected  —  see  text); 
curve  is  extrapolated  to  DOB'S 
for  40 -ton  and  0.5 -ton  events  to 
obtain  comparisons  for  Pre- 
Gondola  II. 
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Fig.  60.  Total  station  energy  (^E)  vs 

yield  after  depth  of  burst  effect 
has  been  removed. 
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Fig.  61a.  Radial  and  transverse  components  at  Station  2BN  for  Pre-Gondola  I  (all 
events). 


Fig.  61b.  Vertical  component  at  Station  2BN  for  Pre-Gondola  I  (all  events). 
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Fig.  61c.  Resultant  velocity  at  Station  2BN  for  Pre- Gondola  I  (all  events). 
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Fig.  62.  Signals  recorded  at  Station  2BN  for  three  different  yields  (0.5-  and  20 -ton 
single  bursts  and  140 -ton  row  charge). 


-55- 


Analysis  of  Results,  Discussion,  and  Interpretation 


TIME  fflSTORY  VARIANCES 

A  comparison  of  the  signals  for  the 
Pre -Gondola  series  of  events  revealed  two 
primary  facts.  The  first  is  that  variations 
in  yield  and  depth  of  burst  had  little  effect 
upon  the  general  appearance  of  the  signals 
at  any  given  station  except  for  amplitude. 
Figures  61  and  62  illustrate  this  fact  very 
well.  These  figures  would  indicate  that 
the  source  varied  only  in  amplitude  and 
very  little  else. 

The  second  fact  revealed  is  that  in 
spite  of  efforts  to  instrument  sites  of 
similar  properties,  the  signals  at  the 
various  stations  bear  little  resemblance  to 
each  other.  For  instance,  the  large  early 
vertical  arrivals  at  Station  2 BN  regularly 
consisted  of  7  or  8  cycles  of  motion  of 
nearly  constant  amplitude  for  all  events. 

At  Station  5N  for  both  the  Bravo  Event 
and  the  140 -ton  row  charge  this  portion 
of  the  record  for  the  vertical  motion  con¬ 
tained  only  a  single  cycle  of  very  high 
amplitude  motion  and  very  little  other 
motion.  Thus  when  peak  velocity  is 
plotted  as  a  function  of  distance  (Figs.  7 
and  8),  Station  5N  consistently  appears 
high.  When  a  least  squares  fit  to  the  data 
is  made,  the  scatter  is  large  and  consid¬ 
erable  uncertainty  exists  about  projecting 
the  results  to  new  sites. 

There  is  no  way  to  tell  whether  the 
differences  in  signals  between  stations  is 
due  to  transmission  path  or  local  effects, 
or  a  combination  of  these. 

Because  of  the  unusual  motion  at 
Station  5N  during  the  Bravo  Event,  when 
PNE-1105^  was  written.  Station  5N 


figures  were  neglected  and  a  high  slope 
for  attenuation  was  given.  When  the  motion 
at  Station  5N  from  the  140 -ton  row  charge 
showed  the  same  peculiarities,  it  became 
obvious  that  new  slopes  should  be  calculated 
for  Bravo.  This  new  value  for  a  (Bravo) 
of  1.54  is  considerably  lower  than  that  re¬ 
ported  in  PNE-1105  (a  =  2.45)  but  shows 
good  agreement  with  the  a  for  the  140 -ton 
row  charge  calculated  for  the  same  stations. 
The  attenuation  constants  of  1.54  and  1.64 
are  lower  than  expected,  but  Station  5N 
strongly  biases  these  numbers. 

As  was  stated,  the  depth  of  burst 
seemed  to  affect  little  except  amplitude 
(Fig.  61).  A  variation  in  yield  for  these 
studies  also  meant  a  variation  in  depth  of 
burst,  thus  an  even  larger  difference  is 
noted  in  the  signals  shown  in  Fig.  62. 

Part  of  the  difference  in  general 
appearance  is  due  to  the  high  level  of 
background  seismic  noise  in  the  records 
from  the  SD-1  Event  which  interferes 
with  the  comparison.  Little  more 
than  these  qualitative  statements  can  be 
made  concerning  these  records  but  they 
are  significant  statements  nevertheless. 

Although  they  are  not  shown,  the 
records  from  the  SD-2  Event  are  almost 
identical  to  those  of  SD-1  except  in 
amplitude.  The  somewhat  surprising  re¬ 
sult  is  that  amplitudes  were  consistently 
higher  for  SD-2,  which  was  decoupled, 
than  for  the  coupled  shot  SD-1,  Upon  re¬ 
flection  this  is  not  so  surprising.  It  is  a 
well  known  fact  that  an  explosively  driven 
metal  plate  will  generate  a  higher  pres¬ 
sure  in  a  target  if  it  is  allowed  to  be 
driven  across  a  gap  of  inches  magnitude 
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as  opposed  to  being  initially  in  contact 

7 

with  the  target.  The  reason  is  that  the 
gap  allows  for  the  potential  energy  of  the 
explosive  to  be  changed  into  the  kinetic 
energy  of  the  metal  plate,  resulting  in  a 
higher  initial  pressure  in  the  target. 

Since  the  explosive  in  the  SD-2  shot  was 
contained  in  a  1/8 -in.  thick  wall  aluminum 
sphere,  4-3/4  in.  smaller  in  radius  than 
the  cavity,  it  is  possible  that  the  result¬ 
ing  higher  pressure  at  the  cavity  wall 
could  cause  a  higher  percentage  of  the 
explosive  energy  to  be  coupled  into  the 
elastic  seismic  signal. 

VELOCITY  IN  THE  FREQUENCY 
DOMAIN 

General  Behavior 

The  CFPV  method  does  only  a  partial 
job  of  spectral  analysis  in  that  it  does  a 
frequency  separation  of  the  signal  but  fails 
to  give  any  information  concerning  the 
total  energy  contained  in  each  frequency. 

The  reason  is  that  a  single  cycle  of  motion 
gives  the  same  value  of  amplitude  as  a 
continuous  sine  wave,  whereas  the  two 
signals  are  vastly  different  in  terms  of 
power,  etc.  Because  of  the  greatly  varying 
time  history  for  different  stations,  the 
CFPV  method  is  unable  to  clarify  to  any 
great  degree  the  relationships  between 
transmitted  signal  and  received  signal. 

A  review  of  Figs.  9  through  25  illustrates 
the  important  capabilities  of  the  CFPV 
method.  A  signal  usually  contains  several 
significant  frequencies  of  motion  related 
to  the  specific  modes  of  wave  propagation 
with  characteristic  wave  periods.  The 
CFPV  gives  an  indication  of  what  these 
frequencies  are,  though  only  one  or  two 
peaks  are  discernible  in  most  of  the  figures. 


In  general  the  vertical  shows  two  peaks; 
one  at  2  or  3  Hz  due  to  late  arrivals,  and 
the  other  at  about  8  Hz  due  to  early  arrivals. 
The  radial  and  transverse  do  not  always 
show  two  main  frequencies.  The  higher 
frequency  is  usually  smaller  in  amplitude 
than  the  lower  and  is  sometimes  missing 
entirely.  The  maximum  motion  is  usually 
found  in  either  the  low-frequency,  late 
arrivals  of  the  radial  wave  or  in  the  high- 
frequency,  early  arrivals  of  the  vertical 
motion. 

The  correlation  of  the  CFPV^s  at  the 
different  stations  for  the  same  event  is 
better  than  for  the  time  history  signals 
themselves.  The  general  shape  of  the 
CFPV's  is  fairly  consistent  except  at 
Station  5N  where  that  single  cycle  pulse 
greatly  distorts  the  CFPV  of  the  vertical. 
This  record  from  Station  5N  illustrates 
one  of  the  main  failings  of  the  CFPV 
method. 

Range  Dependence 

As  was  done  for  the  peak  value  alone,  a 
plot  of  the  peak  value  at  each  frequency  or 
CFPV  vs  distance  (taken  from  Figs.  19 
through  25)  may  be  made  and  an  a  be  cal¬ 
culated  from  a  best  fit  to  the  points.  The 
equation  is 

CFPV(f)  =  V  (f)  = 

(These  plots  are  not  shown.)  When  these 
a’s  are  plotted  against  frequency,  the 
curve  in  Fig.  26  is  produced.  Theory  and 
investigations  have  shown  that  a  should 
be  a  function  that  varies  fairly  smoothly 

g 

and  that  increases  with  frequency.  It  is 
therefore  surprising  to  note  that  the  curves 
in  Fig.  26  are  far  from  smooth  and,  even 
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more  surprising,  they  decrease  with  an 
increase  in  frequency.  This  does  not 
appear  to  be  any  failing  of  the  CFPV  method 
since  the  same  result  is  obtained  by  the 
PSD  method  (see  Range  Dependence  in 
Power  Spectral  Density  section).  It  there¬ 
fore  appears  to  be  inherent  in  the  signal. 

A  possibility  is  that  background  noise, 
both  seismic  and  electronic,  may  be  con¬ 
tributing  factors  but  no  conclusion  can  be 
drawn  until  further  studies  are  made. 

Source  Configuration 

Good  consistency  is  shown  in  the  CFPV 
method  when  the  same  station  measures 
events  at  the  same  SGZ.  This  is  shown  by 
comparing  the  CFPVs  of  the  various 
components  for  Stations  2BN  and  6W  on 
the  SD-1  and  SD-2  Events.  These  are 
shown  in  Figs.  11,  12,  13,  16,  17,  and  18. 

There  appears  to  be  no  significant 
change  in  frequency  peaks  as  the  source 
coupling  is  changed.  But  the  unusual  re¬ 
sult  of  higher  amplitudes  at  nearly  all 
frequencies  for  the  decoupled  event  is 
again  apparent.  It  is  obvious  that  this 
effect  of  increase  in  seismic  signal 
amplitude  occurs  only  for  small  distances 
between  inner  and  outer  spheres  at  the 
sovrce.  Eventually,  as  this  void  space  is 
increased  beyond  some  value,  the  effect 
will  be  to  lower  the  pressure  on  the  cavity 
wall  thus  decreasing  the  seismic  signal. 

Yield  Dependence 

Only  two  points  are  available  to  compare 
CFPV’s  at  different  yields.  The  point  for 
SD-1  (0.5-ton)  is  shown  in  Fig.  9,  and  the 
point  for  the  140 -ton  row  charge  is  shown 
in  Fig.  19,  There  appears  to  be  a  general 
shift  in  the  curves  to  lower  frequencies  at  the 
higher  yields,  though  the  shift  is  small. 


being  only  about  30  to  50%  for  a  factor  of 
280  in  yield.  Such  a  small  shift  may  be 
related  to  the  actual  depth  of  burst,  and 
one  should  be  careful  in  trying  to  extend 
this  to  the  completely  contained  case. 

The  amplitudes  are  different  by  approx¬ 
imately  XIOO.  This  is  only  a  crude  esti¬ 
mate  of  the  ratio  since  an  accurate  one  is 
not  possible  due  to  the  shapes  of  the 
curves.  A  scaling  factor  of  ^  =  0.8  is 
given  by  the  equation: 

Ai/Ag  =(Wq/WP^ 

100  =  (280)^ 

or  A  =  =  kW^ 

where  W  is  the  yield. 

POWER  SPECTRAL  DENSITY  AND 
SEISMIC  ENERGY 

General  Characteristics 

The  PSD’s  for  all  stations  analyzed  are 
much  more  consistent  than  CFPV’s.  This 
is  to  be  expected  since  the  PSD  method  is 
an  averaging  process,  and  the  problems 
exemplified  by  Station  5N  are  overcome. 

Since  it  is  an  averaging  technique,  the 
length  of  record  analyzed  becomes  im¬ 
portant  because  a  single  pulse  will  give  a 
smaller  PSD(f)  for  a  longer  length  of 
record  analyzed.  The  total  energy,  how¬ 
ever,  is  independent  of  the  length  of  record. 
For  this  reason  the  two  following  steps 
were  taken  to  insure  consistency.  The 
total  length  (in  seconds)  of  record  was 
made  the  same  when  possible  so  that 
visual  comparisons  of  PSD’s  would  be  ob¬ 
jectively  correct.  This  length  of  record 
was  32  sec  for  all  stations  and  all  events 
except  for  Stations  2BN  and  6W  on  the 
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SD-1  and  SD-2  Events.  The  latter  records 
were  only  24  sec  long.  Numerical  calcu¬ 
lations  were  made,  whenever  a  choice  was 
possible,  on  integrated  PSD  or  total  energy 
in  a  signal  or  station  (E  or  ^E). 

The  radial  and  transverse  components 
consistently  show  a  pronounced  peak  at 
1  to  3  Hz  with  smaller  peaks  at  6  to  8  Hz. 
The  vertical  also  shows  a  most  pronounced 
peak  at  1  to  3  Hz  and  several  strong  peaks 
at  6  to  10  Hz. 

Source  Configuration 

A  comparison  of  the  PSD’s  for  the  SD-1 
and  SD-2  Events,  component  by  component 
at  Stations  2BN  and  6W,  are  shown  in 
Figs.  31  through  38.  These  show  a  re¬ 
markable  repeatability  for  the  two  events 
with  only  the  amplitudes  varying.  This 
variance  is  again  an  increase  of  signal 
strength  for  a  small  decoupling  factor. 
Enough  has  already  been  said  of  this 
phenomenon  and  the  results  of  the  PSD 
method  are  consistent  with  other  analysis 
techniques. 

Depth  of  Burst 

Figures  47  through  50  show  the  PSD>s 
component  by  component  at  Station  2BN 
for  the  four  20 -ton  events.  These  events 
differed  only  in  depth  of  burst.  Again 
these  figures  show  little  variation  except 
in  amplitude,  a  fact  not  unexpected  since 
the  actual  time  histories  are  so  similar. 
When  the  PSD  is  integrated  over  frequency 
to  obtain  E,  and  E  is  summed  over  com¬ 
ponents  to  give  the  ^E  received  at  the  re¬ 
cording  site,  a  dependence  upon  actual 
depth  of  burst  becomes  apparent.  This 
is  shown  graphically  in  Fig.  59. 

The  datum  point  for  the  Bravo  Event 
shown  in  Fig.  59  was  low,  and  reasons 


for  this  were  discussed  at  length  in  PNE- 
1105.^  Suffice  it  to  say  here  that  there  is 
good  evidence  that  Bravo  vented  early  and 
coupled  less  efficiently  than  it  should  have. 
For  this  reason  it  has  been  neglected  in 
the  following  analysis. 

The  straight  line  in  Fig.  59  is  a  least 
squares  fit  to  the  three  points.  Delta, 

Alfa,  and  Charlie.  It  is  doubtful  that  the 
real  curve  is  a  straight  line,  but  three 
points  do  not  give  enough  clues  to  warrant 
any  other  choice.  As  a  first  approximation 
then,  the  linear  fit  to  a  semilog  plot  is 
made.  For  obvious  reasons  this  curve 
must  flatten  out  at  very  deep  depths  and 
thus  should  not  be  extended  to  containment 
cases. 

It  is  the  contention  of  the  author  that 
this  dependence  of  total  station  seismic 
energy  (^E)  on  depth  of  burst  is  a  result 
of  the  difference  in  time  for  the  expanding 
cavity  to  communicate  with  the  atmosphere. 
This  contact  with  the  atmosphere  causes 
a  sudden  drop  in  cavity  pressure  and,  in 
effect,  an  end  to  its  capability  to  generate 
seismic  signals.  Thus  the  longer  a  posi¬ 
tive  pressure  (i.  e.,  greater  than  over¬ 
burden)  exists  in  the  cavity,  the  longer 
will  be  the  time  during  which  cavity  growth 
and  motion  can  radiate  energy  into  the 
earth  and  hence  the  greater  the  energy 
coupled  into  the  seismic  signals. 

By  extending  the  straight  line  in  Fig.  59 
one  is  able  to  read  directly  the  value  of  E 
at  Station  2BN  for  20  tons  at  any  depth  of 
burst  within  the  limits  of  a  few  meters  to 
nearly  complete  containment. 

Yield  Dependence 

The  variation  of  PSD  with  yield  for  the 
components  and  resultants  at  Station  2BN 
is  shown  graphically  in  Figs.  27  through 
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30.  These  figures  illustrate  the  good 
reproducibility  of  PSD  for  three  different 
events. 

In  order  to  imder stand  the  yield  depen¬ 
dence  of  the  value  of  for  the  cratering 
case,  the  depth  of  burst  effect  must  be 
removed.  This  is  done  by  using  the  re¬ 
sults  of  the  previous  section  and  the  curve 
from  Fig.  59.  The  procedure  is  as  follows. 
Values  of  for  all  events  are  normalized 
to  some  standard  depth  (in  this  case  the 
0.5 -ton  optimum  depth)  by  multiplying  the 
measured  by  the  ratio  of  the  0.5 -ton 
point  to  the  actual  depth  of  burst  point. 

This  new  value  or  ^E»  is  the  ^E  which 
would  result  if  the  event  were  detonated  at 
the  shallower  depth.  The  is  inde¬ 

pendent  of  depth  of  burst  within  the  limits 
mentioned  in  the  previous  section.  If  the 
measurement  were  made  at  Station  2BN 
for  Pre -Gondola  events  this  ^^E’  would 
also  be  independent  of  range.  A  plot  of 
such  events  would  give  the  yield -E 
dependence. 

The  calculated  values  for  the  20 -ton 
events  are: 

=  16.9  X  10“2  (7.7/151) 

=  8.65  X  10"3  (Alfa  Event) 

=  6.34  X  10-2  (7.7/64 

=  7.62  X  10“^  (Charlie  Event) 

=  20.3  X  10-2  (7.7/214) 

=  7.8  X  10-3  (Delta  Event) 

To  calculate  for  the  row  charge  (three 
20 -ton  and  two  40 -ton  charges),  each 
charge  must  be  treated  separately  and 
multiplied  by  its  proportion  of  the  total. 

Thus 
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=  6,64  X  10  ^  (140 -ton  event) 


These  numbers  are  plotted  in  Fig.  60 
along  with  ^^E  for  the  SD-1  Event,  and  a 
Straight  line  fit  to  the  points  is  made  (for 
SD-1,  2^E  =  2!^E>).  The  slope  of  this 
line  is  a  =  1.11.  If  all  of  the  variables 
except  yield  had  truly  been  normalized  out, 
a  slope  of  1.0  or  less  would  be  expected. 
Considering  the  few  numbers  of  points 
used  in  the  depth  of  burst  curve,  this 
close  agreement  is  quite  good.  It  is 
reasonably  safe  to  conclude  that  the  total 
energy  scales  as  the  first  power  of  yield 
in  this  range,  and  that  the  difference  (be¬ 
tween  1.0  and  1.11)  lies  in  the  inaccuracy 
of  the  curve  in  Fig.  59. 

Range  Dependence 

A  comparison  of  the  separate  com¬ 
ponents  at  all  stations  is  shown  in  Figs.  39 
to  42  for  the  Bravo  Event  and  in  Figs,  43 
to  46  for  the  140 -ton  row  charge.  Agree¬ 
ment  of  the  curves  as  a  fvinction  of  range 
is  fair  for  Bravo  and  better  for  the  140- 
ton  event.  The  best  consistency  is  in  the 
resultant  shown  in  Fig.  46. 

A  plot  of  the  peak  energy  density 
(PSD(f)^^^)  for  each  component  and  re¬ 
sultant  vs  range  is  shown  in  Fig.  51  for 
Bravo  and  in  Fig.  52  for  the  140 -ton  row 
charge.  The  lines  are  best  fits  to  the  peak 
resultant  values.  The  attenuation  slopes 
of  4.01  and  3.71  are  in  fair  agreement 
with  each  other. 

An  excellent  agreement  results  between 
the  two  events  when  E  and  ^E  are  plotted 
(Figs.  53  and  54).  The  two  attenuation 
slopes  for  ]^E  are  or  =  3.55  and  q?  =  3.53 
for  the  two  events  (see  the  next  section  for 
a  discussion  of  these  values). 

Attenuation  of  the  energy  in  each  fre¬ 
quency  band  is  not  made  very  clear  by 
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PSD  results.  The  attenuation  slope  for 
each  frequency  is  shown  in  Fig.  55  for 
Bravo  and  in  Fig.  56  for  the  140 -ton  event 
for  three  components.  The  resultants  from 
the  two  events  are  compared  in  Fig.  57. 

Although  there  is  good  agreement  be¬ 
tween  events  it  is  still  not  clear  why  o^Cf) 
is  not  a  simple  function  of  frequency  and 

0 

why  a  decreases  as  frequency  increases. 

Comparison  of  Analysis  Methods 

It  is  clear  from  the  above  discussion 
that  the  CFPV  method  reveals  more  about 
the  nature  of  the  incoming  seismic  signal 
than  a  simple  peak  velocity  determination. 
Likewise  the  PSD  method  reveals  even 
more  information  than  CFPV,  and  the  re¬ 
sults  appear  to  be  more  consistent  with 
the  mean  deviation  of  various  stations  being 
smaller  for  this  latter  method. 

An  interesting  and  encouraging  result 
is  obtained  by  comparing  slopes  to 
peak  velocity  slopes.  Let  us  reason  as 
follows:  (1)  the  peak  velocity,  on  the 
average  (i.  e,,  for  many  stations  and  events) 
will  be  proportional  to  the  signal  attenua- 
tion^  and  (2)  average  energy  is  proportional 


to  average  velocity  squared.  It  then 
follows  that  the  slope  of  the  peak  velocity 
curve  should  be  one  half  the  slope  of  the 
energy  curve.  From  the  foregoing: 

1/2  Of  (  total  energy,  Bravo)  =  1.76 

1/2  a  (total  energy,  140-ton)  =  1.76 
a  (peak  velocity,  Bravo)  =  1.54 
a  (peak  velocity,  140-ton)  =  1.64 
a  (peak  velocity,  average)  =  1.6 

Other  investigators  have  found  a's  for  peak 
velocity  ranging  from  1.5  to  2.0  with  a 
rather  generally  accepted  average  of 
about  1.8.^'®'®  Thus  the  a  of  3,55  for 
energy  attenuation  seems  to  be  an  excel¬ 
lent  number. 

A  comparison  of  1/2  c^(f)  from  the  re¬ 
sultant  PSD(f)  and  a  for  the  CFPV  on  the 
140-ton  event  shows  good  agreement  at 
nearly  all  frequencies  (Fig.  58).  Since 
these  methods  are  different  in  approach, 
the  restilts  support  the  validity  of  each 
method  in  that  they  accurately  describe 
the  content  of  the  signal.  They  do  not  ex¬ 
plain  why  the  high  frequencies  attenuate 
(apparently)  more  slowly  than  the  low  fre¬ 
quencies. 


Conclusions 


From  the  foregoing  the  following  con¬ 
clusions  may  be  made: 

1.  Peak  velocity,  CFPV,  and  PSD  are 
all  valuable  tools  provided  one  under¬ 
stands  their  limitations  and  applicability. 
However,  for  consistency  in  describing 
the  transmitted  signal  in  spite  of  station 
differences,  PSD  is  the  best  of  the  three 
methods.  In  fact,  its  good  consistency  in 
correlating  with  yield  suggests  two 
possibilities: 


a.  This  might  be  a  valuable 
diagnostic  tool  to  determine 
yields  of  undergromd  explosions 

b.  This  technique  might  be  applied 
to  earthquake  records  and  be  the 
basis  for  an  earthquake  magnitude 
scale, 

(In  fact,  I  might  even  be  so  bold  as  to  call 
this  scale  the  "Power  power  earthquake 
magnitude  scale" — albeit,  somewhat 
facetiously). 
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2.  The  total  station  energy  (]^E)  of 
the  transmitted  seismic  signal  over  the 
ranges  of  0.5  to  140  tons  and  probably 
higher  varies  as  the  first  power  of  the 
yield. 

3.  Varying  the  absolute  depth  of  burst 
affects  the  transmitted  seismic  signal  at 
cratering  depths  to  a  marked  extent.  In 
terms  of  energy,  a  factor  of  2  for  a  small 
change  in  depth  of  burst  is  not  unusual  and 
might  be  important  in  selecting  a  depth 
for  a  future  cratering  event. 

4.  The  total  station  energy  (^E)  in  a 
seismic  signal  over  the  intermediate 
range  for  frequencies  of  0.5  to  20  Hz  goes 
as  the  -3.55  power  of  range 


5.  The  effect  of  geologic  differences 
in  transmission  path  and  local  geologic 
columns  causes  gross  differences  in  the 
time  dependent  appearance  of  the  ground 
motion.  These  geologic  differences  are 
not  as  effective  in  altering  PSD  or  the 
total  energy  arriving  at  the  station  when 
the  geologic  and  geophysical  properties  of 
various  stations  are  similar.  A  con¬ 
clusion  as  to  what  would  be  the  result  for 
large  differences  in  transmission  path 
and  local  geologic  properties  cannot  be 
made  on  the  basis  of  this  work,  and 
further  study  is  needed. 
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Appendix  A 


A  Discussion  of  the  Transmission 
of  Seismic  Signals  and  their  Measurement 


The  most  universally  used  method  of 
determining  the  gromd  motion  from  an 
underground  explosion  is  to  measure  with 
a  transducer  a  dynamic  coordinate  of  some 
finite  part  of  the  earth  as  a  function  of 
time.  Ordinarily  the  coordinate  is  dis¬ 
placement,  velocity,  or  acceleration.  In 
theory  any  two  may  be  obtained  from  a 
third  by  integration  or  differentiation.  In 
practice  this  is  difficult  because  of  inherent 
noise,  instrument  drift,  instabilities  of 
numerical  methods,  etc.  It  is  also  usually 
not  possible  to  measure  all  three  directly 
due  to  the  financial  expense  and  instrument 
capabilities.  The  usual  compromise  is  to 
use  the  best  instrument  suited  for  the 
magnitude  of  motions  and  to  do  the  best 
one  can  with  analysis  techniques  available. 
For  the  intermediate  range  (a  few  to  a 
few  hundred  kilometers),  velocity  gauges 
are  the  usual  choice. 

The  time -dependent  motion  that  occurs 
at  any  given  surface  location  is  the  result 
of  a  transmitted  signal  which  has  had  a 
very  complicated,  and  for  the  most  part 
unknown,  past.  This  motion  is  not  one 
signal  but  many,  each  arriving  via  a 
separate  path  from  the  source.  Thus 
there  are  direct  arrivals  traveling  ex¬ 
clusively  in  the  near  surface  material 
while  other  signals  are  reflected  and  re¬ 
fracted  from  the  deeper  layers  of  the 
earth fs  crust.  In  addition  to  various  paths 
there  are  several  modes  of  wave  motion, 
such  as  compressional,  shear,  Rayleigh, 
and  Love  waves  each  traveling  at  different 


velocities.  Finally,  for  each  path,  there 
are  numerous  interfaces  and  impedance 
mismatches  with  attendant  reflections  and 
alterations  in  the  signal.  Certain  geologic 
and  geometric  structures,  particularly 
near  the  surface,  may  act  as  wave  guides 
or  resonating  structures.  Finally,  the 
signal  which  the  investigator  obtains  is 
not  the  ground  motion  but  is  the  response 
of  some  mechanical  system  to  that  ground 
motion. 

The  total  path  a  signal  travels  may  be 
divided  into  four  areas:  (1)  source,  which 
for  our  case  includes  variations  in  the  ex¬ 
plosive  (such  as  yield  and  immediate 
surrounding  area  including  shape  and 
volume  of  source  cavity),  the  immediate 
surrounding  earth  material,  and  the 
relationship  to  the  nearest  interfaces  and 
impedance  changes,  specifically  depth  of 
burst;  (2)  the  so-called  transmission  path 
which  extends  to  all  the  earth  except  those 
areas  immediately  surrounding  the  source 
and  the  recording  stations;  (3)  the  re¬ 
ceiving  site  which  may  include  a  few  to 
several  hundred  feet  of  surface  material 
below  and  around  the  transducer;  and 
(4)  the  transducer  which  includes  itself 
and  the  problems  associated  with  the  re¬ 
cording,  retrieval,  and  reduction  of  data. 

If  one  wishes  to  study  the  effect  of 
some  variable  upon  the  transmitted  signal, 
all  things  should  be  held  constant  except 
that  variable.  This  is  not  entirely  possible 
since  the  so-called  variables  are  not 
completely  independent  of  each  other.  As 
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a  specific  example,  a  sine  wave  source 
of  amplitude  and  frequency  f may 
be  transmitted  through  some  transmission 
path  to  a  receiving  site  with  amplitude 
R^.  If  now  a  new  sine  wave  source  S2  of 
the  same  amplitude  but  new  frequency 
f2  is  transmitted  over  the  same  path 
one  would  expect  a  sine  wave  of  amplitude 
and  frequency  f2  to  be  received  at  the 
distant  point.  This  however,  will  not  be 
the  case  if  transmits  one  frequency, 
say  f^,  better  than  another  or  f2.  The 
transmission  path  may  also  have  amplitude 
dependencies. 

The  problem  is  further  complicated  if 
the  source  is  a  complex  of  frequencies 
and  amplitudes.  It  is  a  little  bit  like  a 
photographer  trying  to  figure  out  how  much 
light  to  put  on  his  subject  when  he  doesn't 
know  what  his  lens  opening  is  or  what 
speed  film  is  in  the  camera. 

It  would  seem  to  be  almost  hopeless  to 
accomplish  anything  against  such  odds. 


However  some  hope  does  exist.  If  the 
nonlinearities  in  the  problem  are  made  as 
small  as  possible  they  might  be  ignored. 
Since  the  Pre -Gondola  series  of  events 
makes  it  readily  possible  to  study  the 
source,  it  is  the  uncertainties  in  the  trans¬ 
mission  path,  local  effects,  and  trans¬ 
ducer  effects  which  are  of  concern. 

The  transducer  is  the  easiest  problem. 
Although  it  is  nonlinear  below  1  Hz,  its 
response  is  well  known  and  understood.  It 
has  already  been  mentioned  in  body  of  the 
report  that  it  was  believed  that  the  geology 
in  northeastern  Montana  is  very  uniform. 

By  choosing  a  radial  path  in  a  direction  of 
little  structural  relief  (in  both  surface 
and  subsurface)  and  by  selecting  station 
locations  over  what  appeared  to  be  similar 
geologic  columns  it  is  possible  to  minimize 
the  effects  of  local  geology  and  trans¬ 
mission  paths. 
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Appendix  B 

A  Discussion  of  Comb  Filter  Peak  Velocity  and 
Vector  Addition  Analysis  Techniques 


A  simple  measure  of  a  seismic  signal 
is  to  measure  the  maximum  amplitude 
without  regard  to  its  frequency.  An  im¬ 
provement  is  to  separate  a  time  history 
into  its  frequency  components  and  to  mea¬ 
sure  the  peak  amplitude  in  each  frequency. 

To  do  this  by  analog  methods  requires 
a  filter  which  passes  only  a  single  fre¬ 
quency  and  eliminates  all  others. 
Mathematically,  the  response  of  such  a 
filter  is  represented  by  a  Dirac-delta 
function.  An  infinite  array  of  these  filters 
would  be  required  to  determine  the  entire 
spectrum  in  a  signal. 

Obviously,  a  device  with  an  infinite 
array  of  filters  does  not  exist,  but  a  finite 
array  can  be  made  to  approximate  it  to  a 
sufficient  degree  of  accuracy.  Filters 
exist  which  pass  a  given  frequency  with  a 
ratio  of  1.0  and,  on  either  side  of  this 
center  frequency,  pass  information  with  a 
ratio  of  less  than  1.0.  A  series  of  such 
filters  (called  comb  filters),  set  at 
appropriate  frequencies,  gives  an  approx¬ 
imation  to  the  actual  spectral  content  of  a 
signal.  This  approximation  is  improved 
as  the  attenuation  on  either  side  of  the 
center  frequencies  increases  and  as  the 
number  of  filters  increases. 

It  is  necessary  to  have  filters  which 
span  the  frequency  content  of  a  seismic 
signal.  Due  to  the  inherent  characteristics 
of  large  explosives  the  meaningful  fre¬ 
quencies  are  below  20  Hz.  Transducer 
characteristics  further  limit  the  recorded 


signals  to  frequencies  greater  than  some 
minimum.  For  the  case  of  a  1 -cycle 
velocity  geophone  the  sensitivity  drops  off 
sharply  below  1.0  Hz  and  seismic  fre¬ 
quencies  of  less  than  0.5  Hz  are  usually 
considered  as  "lost”  to  this  type  of 
instrument. 

Let  us  return  now  to  the  method  of 
measuring  just  the  peak  motion  in  a  signal 
regardless  of  frequency.  A  mechanical- 
electrical  transducer  is  usually  constrained 
to  a  single  degree  of  freedom  in  order  to 
reduce  the  complexity  of  its  behavior. 
However,  a  point  in  space  has  three  degrees 
of  freedom  (three-translational  coordinates) 
and  a  body  in  space  has  six  degrees  of 
freedom  (three  translational  and  three 
rotational  coordinates).  For  motions  in 
the  "intermediate"  range  the  rotational 
motions  may  be  ignored  as  they  are 
infinitesimal  compared  to  the  translational. 

The  standard  transducer  array  measures 
motions  in  three  directions  which  then 
comprise  an  orthogonal,  rectangular 
coordinate  system.  This  means  that  an 
arbitrary  spatial  motion  is  not  measured 
directly.  Each  transducer  measures  only 
the  component  of  motion  parallel  to  its 
degree  of  freedom.  To  regain  the  actual 
motion  one  must  transform  back  to  the 
spherical  coordinate  system  to  obtain  the 
magnitude  (R)  and  the  directional  angles. 
This  is  done  by  the  following  well  known 
relationships 
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system.  Although  this  is  an  exceedingly 
small  error  it  can  be  made  smaller  still 
by  digitizing  x,  y,  z,  y,  xj  x,  y,  z,  y,  x; 
X,  etc.  Thus  the  output  is  x{t),  y(t  +  At), 
z(t  H-  2Ar),  y(t  +  3At),  x(t  =  yAr),  etc. 
Then 

x(t.)  +  x(t.  +  4At)  ^ 

X.  = - 2 -  =  +  2At) 

y(t.  +  At)  +  y{t.  +  3At) 

yi  = - ^ - - 

=  y(t.  +  2At) 


In  all  of  the  above  equations  the  variables 
X,  y,  z.  Res,  0,  and  <})  are  functions  of 
time  and  hence  can  be  operated  upon  jointly 
only  if  time  is  the  same  for  all  variables 
in  a  given  calculation.  This  is  achieved 
readily  if  x,  y,  and  z  are  recorded 
simultaneously  on  the  same  tape.  All  sub¬ 
sequent  operations  are  done  in  a  similar 
manner,  and  the  time  relationship  between 
components  is  held  fixed.  However,  in 
actual  practice  one  complication  arises 
due  to  the  manner  in  which  the  data  are 
digitized  for  these  studies.  The  digitizer 
can  read  only  one  channel  at  a  time;  hence 
X,  y,  and  z  must  be  read  sequentially  as 
the  tape  passes  over  the  head.  This  re¬ 
sults  in  a  small  time  difference  (At) 
between  x^  and  y^  of  0.0005  sec  for  this 


z.  =  z(t  +  2At) 

Since  no  periods  anywhere  near  0.005  are 
present  in  the  data,  x^,  y.,  and  z^  are 
extremely  accurate  (indeed,  much  more 
accurate  than  necessary  in  this  case; 
however,  the  computer  code  used  for  these 
studies  was  written  so  as  to  be  much  more 
generally  applicable).  This  code  was 
written  for  the  IBM  7094  and  provides 
three  outputs: 

1.  A  graphical  plot  of  x.,  y^,  z., 

|Res^|,  (Res)^,  and  cf)^  vs  time 
scaled  to  proper  values. 

2.  A  printed  output  giving  the  actual 
numbers. 

3.  A  digital  tape  suitable  as  input  for 
the  PSD  code. 


-66- 


Appendix  C 

A  Discussion  of  the  Power  Spectral  Density  Code 


/ 


The  power  spectral  density  is  a  mea¬ 
sure  of  the  energy  per  unit  frequency  in 
a  given  time  history  and  is  a  function  of 
frequency.  By  definition  then: 


PSD(f)  = 


AE 

Af 


and  a  measure  of  the  total  energy  is  given 

by 


=  total  energy. 


AT  =  0.01  sec 
At  =  0.0025  sec 
N  =  300 

The  equations  given  by  Blackman  and 
Tukey  are  as  follows: 

The  auto  covariance  is  given  by 
I-i 

\  ~  1  -  *i+je 

i  =  l 

for  i  =  0,  1.  .  . ,  N;  and  I  is  the  total 
number  of  data  points  in  a  record. 

Then  the  power  spectrum  is  given  by 


In  principle  the  PSD(f)  can  be  calculated 
with  integral  methods  (see  Ref.  5  for  a 
complete  discussion  of  theory).  In  practice 
one  can  not  handle  an  integral  directly  by 
digital  methods  since  it  requires  an 
infinite  number  of  points  in  the  calculation. 
There  are  several  techniques  utilizing  a 
finite  number  of  points  which  approximate 
the  PSD,  The  method  used  here  is  the 

5 

one  suggested  by  Blackman  and  Tukey. 

The  accuracy  of  this  method  is  good 
but  care  must  be  taken  in  choosing  the 
lag  length  (AT),  the  number  of  lags  (N), 
and  the  time  between  successive  data 
points  in  a  given  signal  (At)  all  with  re¬ 
spect  to  the  frequencies  present  in  the 
data.  The  data  contain  information 
primarily  between  0.5  and  20.0  Hz,  the 
low  frequencies  being  filtered  out  by  the 
mechanical  response  of  the  system  and 
the  high  frequency  noise  being  electron¬ 
ically  filtered  out.  The  parameters  for 
this  work  done  in  this  study  were: 


N 


+  Ajj  +  2 


N-1 

z 

m=l 


for  k  =  0,  1, .  . . ,  N 


Since  the  above  is  an  estimate,  Blackman 
and  Tukey  recommend  smoothing  by  the 
following  process: 


Pq  =1  (Pq  + 


^  4  ^^k-1  ^k+1^  2  ^k 


Now  P^  is  the  approximate  energy  in  a 
given  frequency  band.  To  obtain  the 
energy  density  or  PSD(f)  we  must  divide 
by  the  Af  or  frequency  band  width, 
where 
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Then 


PSD(f)  =  (2ATN) 


The  total  energy  in  a  component  then  is 


E.l  PSD(Odt»^ 


IK 


A  measure  of  the  total  energy  for  all 
three  components  is  just  the  sum  of  E  or 


for  i  =  X,  y,  and  z. 


or 


Since  the  PSD  is  calculated  on  three 
mutually  orthogonal  coordinates  (x,  y, 
and  z)  the  energy  in  each  coordinate  is 
mutually  isolated  and  independent  of  the 
other  two  coordinates.  Thus  one  method 
of  calculating  the  resultant  PSD(f)  is 
simply  to  sum  three  components  or 

for  i  =  X,  y,  and  z. 


N 

■  Z  At. 

k=0 

The  code  (PSD)  developed  to  do  the 
above  calculations  was  written  for  the 
CDC  3600  computer  and  utilized  a  tape 
written  by  the  IBM  7094  from  the  VECTOR 
code.  The  output  is  in  two  forms:  visual 
(PSD(f)  vs  frequency  plots)  and  digital 
(numeric  printout). 
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